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ABSTRACT 


A search was made for 183 GHz line emission from water vapour m the 
direction of twelve Mira and two semiregular variables. Upper limits to the 
emission are in the range of 2000 to 5000 Jy. It is estimated that thermal 
emission from the inner regions of late-type stellar envelopes will be on 
the order of ten Jy. Maser emission, according to one model, would be an 
order of magnitude stronger. From the limited set sampled, the possibility 
of very strong maser emission at 183 GHz cannot yet be ruled out. 



I. INTRODUCTION 


Molecular line spectra at centimeter and millimeter wavelengths have 
proven useful in determining the conditions in the envelopes surrounding 
giant and supergiant stars of late spectral type. They have shed light on 
the gas-phase composition, isotopic abundances, radial structure of 
temperature and density, and mass loss rates (e.g., Zuckerman 1980a, b). 

While maser transitions are more difficult to interpret than transitions 
between thermally populated levels, their sensitivity to particular 
combinations of environmental parameters gives them the potential of being 
rather fine probes of specific regions around a star. Water vapour emission 
at 22 GHz arises in the acceleration region of the envelope within a few 
stellar radii of the star (~ 10 1 ** cm), where the temperature is 1000-2000 K 
(Vardya 1970; Goldreich 1980). SiO masers also occur in the inner region of 
the envelope. 1612 MHz OH masers occur at > 1000 stellar radii (> 10^ 6 cm; 
e.g., Bowers, Johnston, and Spencer 1981), where the temperature is ~ 100 K 
(Elitzur, Goldreich, and Scoville 1976). 

Another potential maser probe of stellar envelopes is the 183 GHz line 
of water. In modelling the excitation of water in circumstellar envelopes, 
Deguchi (1977) has predicted an inversion of the 183 GHz transition. That 
this transition does indeed invert under suitable conditions is demonstrated 
by the intensity of the wide wings of this line seen toward Ori KL (Kuiper, 
Zuckerman, and Rodriguez Kuiper 1981), and the variability of the emission 
from this source (Kuiper et_ al_ 1984). However, the conditions in the region 
in Orion where this emission arises are believed to be quite different from 
the stellar conditions modelled by Deguchi. A credible a priori prediction 
of 183 GHz emission in any actual region is probably out of the question 
because the physical parameters such regions are not sufficiently well known 
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(Goldreich 1980). We were therefore motivated to make a search for possible 
183 GHz maser emission in the envelopes of late-type giants and supergiants. 

Because of the limited amount of observing time available at our 
observing altitude, we chose a strategy of surveying a large number of stars 
for strong emission, spending between one half and one hour on each star. 

We selected the stars primarily on the basis of OH and 22 GHz H 2 O maser 
emission. Eleven of the fourteen stars we observed have such maser emission 
(Klemmann, Dickinson, and Sargent 1978). This criterion selects late M 
stars with [0] > [C] (Vardya 1970; Wyckoff and Clegg 1978), and therefore 
those with generally lower mass-loss rates (Zuckerman 1980; Knapp et_ al_ 
1982). A number of preferred candidates could not be observed because they 
were not visible above 35 deg elevation at night on our assigned flight 
dates, because they could not be included in an efficient flight plan, or 
because of equipment failure. The stars observed are listed in Table 1. 

II. OBSERVATIONS 

The observations were made m 1980 February and August with a cooled 
Schottky-diode receiver mounted on the 91-cm telescope of the G. P. Kuiper 
Airborne Observatory. The equipment used and the details of the calibration 
are extensively described in Kuiper et_ £l_ (1984). We present here only a 
summary and additional facts of particular relevance to the observations 
described in this paper. 

The system temperature was ~ 900 K, double sideband. The aperture 
efficiency was 0.40, giving a sensitivity of 10** Jy/K. Typically, an r.m.s. 
noise level of 0.1 K was achieved in half an hour. Pointing was done by 
tracking the star observed, or an adjacent star, in a co-aligned tracking 
telescope, and is estimated to be accurate to ~ 10 arcsec, which is a 
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negligible fraction of the 8 arcmin beamwidth. 

Two spectrometers were used. A 36-channel filter bank has 1-MHz 
channels covering the central 11 MHz, twelve 3-MHz channels covering the 
ranges from 4.5 to 22.5 MHz above and below the receiver’s center frequency, 
and ten 8-MHz filters covering the ranges from 12 to 52 MHz relative to the 
center. (In addition, there were three broadband channels.) The other 
spectrometer was a digital FFT device, covering the central 10 MHz with 256 
channels, for any narrow spectral features that might be detected. 

The principal obstacle to achieving the theoretical sensitivity of the 
receiver for broad lines proved to be severe standing waves in the telescope. 
To minimize these as much as possible, we observed by chopping the 
subreflector in azimuth, symmetrically about the optical axis of the 
telescope. The total throw of the telescope beam was two beamwidths . The 
optical axis of the telescope was positioned one beamwidth either left or 
right of the star being observed. The signal received from the direction 
of the star and that from the blank sky on the opposite position of the 
optical axis were subtracted synchronously in the spectrometers. In the 
difference spectra, the amplitude of the standing wave was ~ 10 K in 1980 
February. The dominant periodicity was ~ 70 MHz, which corresponds to a 
reflection path of 7 m, twice the distance between the receiver and the 
subreflector. In addition, the pattern often appears to be partially 
modulated by a second sinusoid with half-period of ~ 100 MHz, which 
corresponds to the distance between the receiver and a bulkhead port in 
the telescope cavity through which the beam from the receiver must pass to 
reach the tertiary flat. While the standing wave pattern consisted mainly 
of a dual sinusoid, the pattern was not regular enough so that it could be 
entirely removed by fitting a function. However, the pattern was smooth 
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enough that our sensitivity to narrow (~ 1 MHz) lines was not significantly 
degraded. In 1980 August, we reduced the dominant standing wave pattern to 
about half by gluing to the center of the subreflector a small disk of 
Eccosorb foam with a diameter equal to the subreflector's image on itself. 

The results are presented m Table 1. The first, second, and third 
columns give the name, spectral type, and variability type of each star 
observed, taken from Becvar (1959), Gehrz and Woolf (1971), and Kukarkm 
(1969). The fourth column contains the LSR velocity at the center of the 
spectrometers. The effective velocity coverage, with a resolution of 3 MHz 
(4.8 km s - - 1 ), was 45 MHz (73 km s -1 ) centered on the LSR velocity to which 
the receiver was tuned. The upper limit to the 183 GHz flux, given in the 
fifth column, refers to narrow spectral features (~ 1 km s _1 ) which we 
consider the more likely to be detected (see Discussion). It corresponds to 
twice the theoretical r.m.s. noise level. From our flat spectra, such as 
that of V1057 Cyg in Figure lb, we were able to verify that we did obtain 
a channel-to-channel noise level consistent with our system temperature. 
Because of the uncertainty in our baselines, we consider this also to be 
the detection limit for broad spectral features. 

We feel that the 2-o upper limit also applies to narrow features in 
spectra having poor baselines. Figure la shows a spectrum of a Her in which 
such a feature appears to be present. The feature also appeared in the high 
resolution digital spectrometer. We cannot be confident that this is a real 
signal because the observation of this source was hindered by mechanical and 
tracking problems with the telescope. The spectrum consists of only one 
left-right pair, of which both parts were terminated early due to loss of 
tracking. The spectrum has the shortest integration time and the most 
irregular baseline of all our observations. The amplitude of the feature 
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is only 2-3 a, depending on where one chooses to draw the baseline. It does 
illustrate, however, that narrow features could be recognized even m the 
presence of strong standing waves. No other narrow features as prominent as 
the one m the a Her spectrum were observed. 

In the spectra with poor baselines, our ability to detect broad spectral 
features is reduced. Such spectra have been identified with a cross next 
to the flux upper limit in Table 1. 


III. DISCUSSION 

In order to assess the significance of the results, we can estimate 
the intensity of the signals that might be expected. 

We consider first possible thermal emission. The flux density from 
water in local thermodynamic equilibrium m a stellar envelope at a distance 
D, in a line which has a width of AV, may be written as 


x R 2 0 N env f 3( 1 , 3) Ahv 
4TTD 2 (vAV/c) 


(la) 


where Xy g is the relative abundance of water; N gnv is the total number of 
gas molecules in the envelope, * 6 x 10 5 times the mass of the envelope in 
£3(1,3) is the number of molecules in the 31,3 rotational state of 
water, and is given in Eq. 2 of Waters et_ ad . (1980); A is the probability 
of spontaneous emission, * 3.6 x 10 -6 sec -1 ; and h, v, and c have their usual 
meanings. The equation may be conveniently expressed as 



7.4 x 10 11 X, 


h 2 o 


*tnv 


D 

-2 

AV 

KJ 


kpc 


km s -1 


T-l-5 exp (- — ) 
T 


(lb) 


The masses of late-type stellar envelopes have been estimated in the range 
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of a few tenths to a few solar masses. The abundance of water is predicted 
to be almost as high as 10~ 3 because all the available oxygen not bound to 
carbon is m the form of water (Goldreich and Scoville 1976). However, we 
could imagine it to be as low’ as 10” 6 if interstellar cloud chemistry is 
typical (e.g., Prasad and Huntress 1980). Adopting an appropriate average 
temperature presents some difficulty since the temperature in such envelopes 
varies from ~ 10 3 K near the star to a few tens of K, or less, near the 
outer part of the envelope. The radial dependence is probably steeper than 
inverse square root in the inner envelope because of absorption by dust in 
the envelope. We will adopt a temperature of 100 K because the emissivity 
of the 183 GHz line is not a strong function of temperature, until the 
temperature drops below 60 K (see Fig. 2 in Waters et_ al_ 1980). The flux 
density from 183 GHz line emission can then be expressed as 


S(LTE,max) *= 


3200 


' M env' 


O 

1 

CsJ 

X(H 2 0) 


. 


Q 

i 


io" 3 



30 km s~l 


AV 


Jy. 


(2a) 


Assuming the lifetime of the outflow to be ~ 10 4 y (Knapp et_ al_ 1982), we 
can estimate the mass of the envelope from the mass-loss rate, and write 
the equation as 


S(LTE,max) 


320 


dM/dt 


kpc 

2 

x(h 2 o) 


30 km s' 1 

10" 5 l^yr" 1 


D 


10" 3 


AV 


Jy. 


(2b) 


This level of thermal emission, however, cannot be expected when one 
considers the density requirement for collisional excitation. Collision 
rates have been calculated by Green (1980) to be in the range of 10”* 2 to 
10"“ cm 3 sec _i for the 3^ and 2 2 q levels over a wide range of temperatures. 
The probabilities of spontaneous emission via the most probable pathways are 
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0.13 sec -1 and 0.27 sec -1 , respectively, for these levels. Thus, if the 
stellar envelopes were optically thin in these transitions, densities on the 
order of 10*° cm -3 would be required for collisional excitation to maintain 
the level populations against spontaneous decay. The density requirement is 
reduced if the emitted photons are reabsorbed by the envelope, since only 
photons which escape from the envelope cause a net de-excitation. This 
trapping cannot be very significant, however, since the expansion of the 
envelope causes most potentially absorbing molecules to be Doppler shifted 
away from the rest velocity in the frame of the emitting molecules. The 
radial density dependence of a uniform mass-loss may be expressed as 


n(H 2 ) - 1.5 x 10 5 


dM/dt 


Thus, we see that collisional excitation can be expected to occur only 
within a few stellar radii (a few times lO* 4 cm) of the surface for a 
mass-loss of 10~ 5 yr *. Since, in a uniform outflow, the mass fraction 
interior to a given radius depends linearly on the radius, and since 
envelopes may extend to a few x 10 i!s cm, we can only expect the few inner 
percent, or less, of the envelope to emit thermal 183 GHz emission. The 
intensity of the line from the region is also reduced by the temperature 
dependent factor (see Eq. lb) because the inner core of the envelope has a 


temperature of ~ 1000 K. In the sixth column of Table 1, we have listed 


the distances of the stars, taken from Bowers, Johnston and Spencer (1983), 
Humphreys (1975), Knapp et_ al_ (1982), Morris and Jura (1983), and Wilson 
et al (1972). In the seventh column we have tabulated estimates of the 


mass-loss rate, taken from Bowers, Johnston and Spencer (1981), Gehrz and 
Woolf (1971), Knapp et_ al^ (1982), and Morris and Jura (1983). From these 
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we have calculated, and tabulated m the eighth column, the factor 
(dM/dt)D~ 2 in convenient units to facilitate comparison of the measured 
flux limits with Equation 2b- To allow for the expected line width m Eq. 
2b, we have also tabulated, m the ninth column, the outflow velocity, 
taken from Bowers, Johnston and Spencer (1981), and Knapp e£ al_ (1982). It 
is clear that the present observations are about three orders of magnitude 
too insensitive to detect thermal 183 GHz water emission from these stellar 
envelopes - 

Deguchi (1977) has modelled the excitation of water m the inner 
envelope of a prototypical mass-loss star, in the range between one and two 
stellar radii (assumed to be 6 x 10 i3 cm). The velocity was assumed to be 
proportional to the distance from the surface of the star, the temperature 
dependence inverse square-root, and the mass-loss rate to be 1(T 5 M yr -1 . 
He computed a 183 GHz maser power of 1.3 x 10 43 photons s _i . Scaling his 
result to the mass-loss rate, it corresponds to a flux at the earth of 


s 183 


dM/dt 


kpc 

2 

km s -i 

io* 5 ^yr _i _ 


D 


AV 


(4) 


In Deguchi’ s velocity model, the amplification occurs in the acceleration 
region so that velocity coherence occurs only where the flow is perpendicular 
to the line-of-sight . Thus, the power may be expected to be confined to a 
thermal linewidth of 1 MHz (■ 1.6 km s -1 at 1000 K) centered at the LSR 
velocity of the star. For a star at 0.2 kpc with a mass-loss rate of 10~ 5 
yr -i , the expected flux would be 360 Jy. Thus, the flux from Deguchi's model 
star at 200 pc distance would lie about an order of magnitude below our 
detection level. Scaling from Deguchi's model, the stars observed should 
have emission about two orders of magnitude weaker than our detection level. 
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The possibility that some stars may have intense maser emission well 
m excess of Deguchi's model, and potentially detectable by us, cannot be 
ruled out on the basis of our limited survey which sampled only 13% of the 
known 22 GHz water maser stars (Kleinmann, Dickinson, and Sargent 1978). 

Prospects for future searches should be quite good. With current 
technology, a Scottky-diode system with 400 K SSB system temperature is 
possible at 183 GHz. An SIS receiver might achieve half that, giving an 
order of magnitude better sensitivity than the search reported here. 

Careful attention will need to be paid to reducing standing waves m the 
telescope system. Path-length modulator (e.g., a reciprocating mirror; 
Gustincic 1977) has been found to be effective in cancelling the effects of 
reflected radiation (Payne, 1983). Another order of magnitude improvement 
in sensitivity can be achieved by using a 3— m balloon-borne telescope 
(Hoffmann, Fazio, and Harper 1983; Melchiorri 1983). A combination of these 
two advances will permit us to detect or put useful limits on 183 GHz maser 
emission from late-type stars. Eventually, large orbiting telescopes such 
as FIRST (Olthof 1983) and LDR (Swanson et_ al_ 1983) should permit the 
detection of thermal 183 GHz line emission from water in stellar envelopes. 

We express our gratitude to P. Batelaan and W. Ricketts for the 
engineering support which made this work possible, to M. Klein for assistance 
with the February observations, and to the staff of the G. P. Kuiper Airborne 
Observatory for their unstinting technical and observational support. The 
research described in this paper was performed at the Jet Propulsion 
Laboratory, California Institute of Technology, under contract with the 
National Aeronautics and Space Administration. 
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Table 1. Results of Search for 183 GHz Line Emission 


Receiver 



Spectral 

Var. 

Central 

max 



(dM/dt )/D 2 



Star 

Type 

Type 

V lsr 

s 

D 

dM/dt 

Vo 

Refs 




(km 8 _1 ) 

(Jy) 

(kpc) 

% y r_1) 

(kpc 2 )/(10~ 5 

% y r_1) 

(km s _l ) 


R Aql 

M5e-M8e 

M 

50 

2000 

0.3 

8xl0 -7 

0.89 


4 

RR Aql 

M6e-M7e 

M 

31.5 

2400 

0.4 

4xl0~ 7 

0.25 


2,6 

R Cas 

M7e 

M 

20 

1300 

0.2 

6.6xl0~ 7 

1.65 

11 

8,11 

U CVn 

M7e 

M 

21 

3300 a 






S CrB 

M6e-M8e 

M 

1.5 

3300 a 

0.4 

4xl0 -6 

2.5 


4,9,11 

NML Cyg 

M7I 


0.5 

2000 a 

2.0 

6.4xl0 -5 

1.6 

23 

2,7,10 

V1057 Cyg 



4 

2000 


9xl0 _e 




a Her 

M51I 

SRc 

0 

5000 a 




4 

U Her 

M6.5e-M8e 

M 

-14 

4000 a 

0.3 

2.6xl0 -6 

2.89 


4,11 

R Leo 

M6 . 5e-M9e 

M 

0 

2700 a 

0.3 

8.5xl0~ 7 

0.94 

7 

8,9 

R LMi 

M7e 

H 

1.5 

2700 

0.4 

10~ 6 

0.63 

6 

8 

R Lyn 

S3e 

M 

2 

4000 a 





1 

R Peg 

M7e 

M 

26 

2000 





1 

VX Sgr 

M41e-M9a 

SRb 

0 

3000 a 

1.5 




3,5,9 

a Sensitivity limited by 

standing waves . 
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FIGURE CAPTIONS 


Figure 1. (a) The spectrum observed in the direction of a Her in the 183 

GHz line of water vapour. The central channels are 1 MHz wide 
those adjacent, 3 MHz; the outer channels, 8 MHz. The full 
length of the error bar on the left corresponds to twice the 
r.m.s. noise level for the 1 MHz channels as calculated from 
the measured system temperature, (b) As above for 1057 Cyg. 
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